Permanent magnetic moment in mesoscopic metals with spin-orbit interaction 
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We argue that at zero temperature an isolated metal particle (or an AB ring) with spin-orbit 
interaction and odd number of electrons will have a permanent magnetic moment, even in zero 
magnetic field (fiux). In a zero-field-cooled state both the direction and the magnitude of the 
' moment varies from particle to particle and averages to zero. In a field-cooled state it averages to 

I ~ Ats (kpl)^^^- We argue that the permanent moment is due to an uncompensated electron in the 

^SJ ■ last occupied (Fermi) level. We introduce an effective single-electron Hamiltonian which accounts 

for spin-orbit coupling. 

OO : I. INTRODUCTION 

Permanent currents in mesoscopic SNS junctions, that is currents flowing in the absence of the phase difference 
between the superconductors, were predicted by Altshuler and Spivak ||l[ in the circumstance of a time-reversal 
^2 l symmetry breaking perturbation, such as spin-flip scattering. Since their model is equivalent to that of an Aharonov- 
Bohm (AB) ring threaded by a flux, it also implies a permanent AB (persistent) current (or, in general, a permanent 
do ' magnetic moment) when the time-reversal symmetry is broken. Such currents (moments), however, average out to 
^ . zero and are uncorrelated with the flux (fleld). 

^ ' Kravtsov and Zirnbauer (KZ) Q predicted permanent AB currents in rings with odd numbers of electrons and 
(non time-reversal breaking) spin-orbit (SO) scattering. Their argument was based on the Kramers' theorem in 
the context of disordered mesoscopic systems Namely, the level degeneracy should remain two- fold, where a state 
and a time-reversed state are two different states at the same energy. KZ realized that in the presence of SO coupling 
these states are current-carrying. To estimate the current's magnitude for strong SO they used the level correlation 
function of a symplectic ensemble 1^,(1). The letter contains a peaked function of the flux Q , Q which, in the case 
Q of a zero flux, reduces to a (5— function. The physical interpretation suggests that the double degeneracy is lifted by 
^ the flux and one of the two current-carrying states becomes preferable; in zero flux either state can be occupied with 
, equal probability. 

^ ' In the absence of SO scattering the two-fold Kramers' degeneracy for an odd number of electrons is trivial in that 
^ ] it corresponds to the spin-up and spin-down states of an uncompensated electron in the last occupied (Fermi) level 
(all other levels being doubly occupied). For a sufficiently strong SO, spin is no longer a good quantum number and 
^\ one must consider an effective moment, whose origin is in mixing of orbital and spin degrees of freedom and which is 
conserved |3| since a metal particle or a ring can be viewed as a closed system. However, unlike atomic physics, for 
instance, where the nature of such moment is obvious as J = L + S, in a disordered conductor L is not a good 
quantum number. It is clear, on the other hand, that the moment is due to the last uncompensated electron, as in 
spin magnetism. 

Since one must differentiate between the odd and even numbers of electrons, understanding of the nature of 
' permanent moments requires formalism that accounts for a single-electron contribution. The purpose of this work is 
' to introduce a single-electron Hamiltonian that would effectively incorporate SO coupling. 
I In a previous paper we addressed a more rigorous derivation of the permanent AB current in the SO case. 

Towards that end we applied a so called "mixed" approach where the Imry's formula jl^ 
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5t ' relating the free energies of the canonical and grand canonical ensemblesfj is used in combination with the exact 
formulae H, ||] for the disorder-averaged level-density correlation function {6viSv2)- We call this approach "mixed" 
because (|l|) implies large particle number fluctuations in the equivalent grand canonical ensemble, (i57V^) ^ 1, and 



*serota@phy sics.uc.edu 

'^In a closed Kubo particle the total number of electrons is presumed fixed, hence one must derive the free energy of a canonical 
ensemble 
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consequently T ^ A, where A is the mean energy-level spacing. In such a circumstance a perturbative expression 
for the level correlation function should, in principle, be used. However, as will be discussed below, the "mixed" 
approach gives qualitatively correct results on approaching the level-quantized limit for T ~ A, when the single- 
electron effects become important. 

In what follows, we will first compare the results obtained from in the "mixed" approach for spin and orbital 
magnetism with those obtained in the perturbative limit and in a single-electron approximation in the level quantized 
limit §, H- [||. Such comparison will enable us to gain insight into the level-quantized problem in the SO case by 
inference from the "mixed" approximation. We proceed to introduce a single-electron Hamiltonian which effectively 
includes SO coupling. We discuss the physical consequences of our ansatz and check them against other known results, 
such as mesoscopic fluctuations. 



II. PERTURBATIVE, "MIXED" AND SINGLE-ELECTRON APPROXIMATIONS 



The results in this Section are derived in our previous papers @, IHl- For simplicity, we consider only 2D 

systems. We begin with the spin magnetic susceptibility and use the notations 



XP = Ms (w) and xc 
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for the Pauli and Curie susceptibility respectively, where {v) is the mean density of levels at the Fermi surface and 
{v) A = s = 2 is the level degeneracy. The summary for the mean susceptibility is given in Table I: 

TABLE I: Spin magnetic susceptibility 



Pertubative 


"Mixed" 


Single-electron 


r> A 


T - A 


T< A 


XP + ^4 [t) J 




const 


XC 

xp{^r 


odd number of electrons 
even number of electrons 



(4) 



Clearly, as the temperature approaches the level spacing, the "mixed" approach correctly predicts the onset of the 
Curie-like susceptibility. Conversely, it is beyond its range of applicability to distinguish between the odd and even 
numbers of electrons and give a single-electron contribution. (In a sense, it averages between the odd and even 
contributions which is dominated by the odd-case Curie susceptibility.) For further comparison with the SO case, we 
note that the "mixed" result is due to Zeeman splitting of the ^— function in the exact level-correlation function. We 
also note that the linear response regime for the spin susceptibility extends to the fields hbH ^ A, at which point the 
energy levels begin to separate into two non-interacting series of levels with spin-up and spin-down electrons Q and 
the magnetic component of the free energy (that is, the change in free energy due to Zeeman splitting) approaches 
-A. 

(We briefly discuss the case of even number of electrons 1 16|. The constant and the exponent a in Table I depend on 
the symmetry a = I for Gaussian orthogonal ensemblefCOE) and a = 2 for Gaussian unitary ensemble (CUE). 
Notice that CUE was not considered in [|l^ since it was assumed that the transition from GOE to GUE takes place 
at plbH ^ A when, in reahty, it takes place at ^ibH ~ Ay/ A/Ec <€. A (see eq. (^) below and Refs. jj]). On the 
other hand, for a strong SO, the spin susceptibility in the Gaussian symplectic ensemble (GSE), considered in | p6[ , is 
not particularly meaningful since spin is not a good quantum number in this case). 

We next turn to orbital magnetism in the absence of SO interaction. This case is summarized in Table II in terms 
of the mean AB (persistent) current !{(/)): 

TABLE II: Persistent current in the absence of SO 



Pertubative 


"Mixed" 


Single-electron 


r> A 


A 


T< A 


gC/^c T 4>o 271-3 CiA ^ 


SeE,^ ^eA^ 

^ (pQ Ztt-^ <p 

(t><.(t>c <l>c <- 4> <- 4>o 


- 2TTeEc4- 
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where 4>o — 21x1 e is the flux quantum (in the units where h - 
coefficient and L is the circumference of the ring. The flux scale ^ 



c — X) and i?c — Djl? where D is the diffusion 
c defines the limit for linear response and, at T = 0, 
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sets the scale of the transition from GOE to GUE ||]. We point out that the "mixed" approach correctly predicts 
current saturation as T crosses A, which is confirmed by a single-electron type calculation where the orbital 
magnetic susceptibility is attributed to the van Vleck response involving virtual transitions between the last occupied 
and the first unoccupied states at the Fermi level. We also note that, as crosses (pc the maximal current is 
^ eA(fc_F£) ' (corresponding to a magnetic moment of order \iB (-E'c/A) ' ~ [Ib (kr^) ), where kp is the Fermi 
wave vector and £ is the electron mean-free-path, and the magnetic energy —I{(t>) 4> saturates to —A. 
We finally turn to the SO case. The results are summarized in Table III: 



TABLE III: Persistent current with strong SO 



Pertubative 


"Mixed" 


Single-electron 


T> A 


A 


r< A 


IpE AA ^pAiSL 

< 0c 0c < 0< 00 


< 0c 0c < < 00 


I eE ^-^ 



(7) 



where the question mark indicates the result inferred from the "mixed" approach. A formal extension to T = in the 
mixed approach yields a permanent current for = Q , 

/(0) = ^v/2^ (8) 

and the same order of magnitude for a permanent AB current is expected in the single-electron picture (and a 

corresponding magnetic moment of order (Ec/A)^^^ ^ fiB [kpf-)^^'^)- In the next section we introduce a single- 
electron Hamiltonian which indeed reproduces this result. 



III. EFFECTIVE SINGLE-ELECTRON HAMILTONIAN 



A single-electron Hamiltonian which effectively accounts for multiple SO scattering effects can be written as 

He// = AL-S + A*BH-(L + 2S) (9) 



where 




[kpf-Y 



(10) 



In the 2D case considered here and magnetic field perpendicular to the plane of the sample, one need to consider only 
z— components of the spin and angular momentum. Treated as a perturbation to second order, the Hamiltonian (|^) 
yields the energy shift given by 

,5^; - -/iB2 (1 + AA,,) S,H - /i|A,,i/2 _ x^^,,Sl (11) 

where 



(12) 



The first term in ( [ll|) describes the coupling of the permanent moment - a combination of spin and SO-induced 
components - to the magnetic field. The second term is the van Vleck term due strictly to orbital effects. The last 
term describes the single-axis anisotropy. In eq. (|l^) we limit our consideration to the virtual transitions between the 
last occupied (Fermi) level i and the first unoccupied level /, so that Kzz > 0. 



In disordered conductors the orbital magnetic effects are large due to diffusion 1 18|. In the present case, it is reflected 
in largeness of A which was evaluated in ||l^ 
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MS [kp^f''] + const2 [\xl\ (kpi) H] (14) 



Consequently for the magnetic moment we find 

dSE 

M = — - — = consti 
dH 

where xl is the Landau susceptibihty | p2[ . The first term in ( p^ is the permanent magnetic moment due to SO 
coupling and the second term is the induced van Vleck orbital moment. 

Independent checks of ansatz (|9|) come from the consideration of the anisotropy term (the third term in (pi])) and 
the mesoscopic fluctuations of spin-polarization ||l^. For the former, using eqs. (|l0| ) and (p^, we find 5E ~ —A, 
which is expected for the energy associated with the transition to GSE. For the latter, the large fluctuations of local 
magnetization were interpreted in terms of the enhancement of the effective electron 17— factor, 17 = 2 (1 + AA^z), to 
1/2 

~ {kpt) , which is consistent with the present results. Another consistency check comes from the requirement that 
both the linear and quadratic magnetic terms in ([l^) be less than A and we find that either one is equivalent to the 
condition (f) < (pc for the transition to GUE, where (pc is given by eq. (|^). 

IV. DISCUSSION 

We proposed the single-electron Hamiltonian (^) which effectively incorporates the effect of multiple SO scatterings 



for the electron in the Fermi level. The resulting magnitude ( 14 ) of the permanent magnetic moment of a 2D ring or 

a simply-connected conductor is consistent with the results based on the perturbative and "mixed" approaches. We 

1/2 

also reproduce the correct energy scales and the enhanced magnitude ~ {kpt) ' of the electron factor. However, 
at present there is no microscopic derivation of (^). We hope to address this in a future work. 
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Note added: After the completion of this paper I became aware of the manuscript ;ond-mat/0001431 which discusses 
some of the same issues. For 2D systems there is a general agreement between the two and with our previous work 
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